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Abstract 

To gain insight into the physics of the metal insulator transition and the effectiveness of cluster 
dynamical mean field theory (DMFT) we have used one, two and four site dynamical mean field 
theory to solve a polaron model of electrons coupled to a classical phonon field. The cluster 
size dependence of the metal to polaronic insulator phase boundary is determined along with 
electron spectral functions and cluster correlation functions. Pronounced cluster size effects start 
to occur in the intermediate coupling region in which the cluster calculation leads to a gap and 
the single-site approximation does not. Differences (in particular a sharper band edge) persist in 
the strong coupling regime. A partial density of states is defined encoding a generalized nesting 
property of the band structure; variations in this density of states account for differences between 
the dynamical cluster approximation and the cellular-DMFT implementations of cluster DMFT, 
and for differences in behavior between the single band models appropriate for cuprates and the 
multiband models appropriate for manganites. A pole or strong resonance in the self energy is 
associated with insulating states; the momentum dependence of the pole is found to distinguish 
between Slater-like and Mott-like mechanisms for metal insulator transition. Implications for the 
theoretical treatment of doped manganites are discussed. 

PACS numbers: 71.10-w,71.30.+h,71.38.-k 
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I. INTRODUCTION 



Correlation-driven metal-insulator transitions are one of the fundamental issues in elec- 
tronic condensed matter physics [1] . It has known for many years that in models of electrons 
moving on a periodic lattice and subject to short-ranged interactions, insulating phases may 
occur if the number of electrons per unit cell, A^, is an integer or rational fraction. If 
is an even integer, then the Pauli principle makes insulating behavior possible even in the 
absence of electron-electron interactions (this is the case of filled bands). If is not an 
even integer, interactions are required to obtain an insulator. Interactions may lead to a 
spontaneous breaking of translational symmetry, implying a unit cell larger by a factor of A 
than the crystallographic cell. If A is such that AA^ is an even integer, insulating behavior 
is clearly possible; this is referred to as the "Slater mechanism" On the other hand in- 
sulating behavior may occur for A^ not an even integer, even in the absence of spontaneous 
translational symmetry breaking. This is often termed the "Mott mechanism" 3|. 

Important new insights into metal-insulator transition physics have come with with the 
development of dynamical mean field theory (DMFT) Single-site dynamical mean field 
theory, which entirely neglects spatial correlations in the electron gas, is found to predict 
a zero-temperature metal insulator transition as the local interaction U exceeds a critical 
strength f/««"9'e-site density A^ = 1 41. More recently, cluster dynamical mean field 

n n Q 

theories [5|, |6|, |7|] which include some degree of intersite correlations have been shown to 
produce metal insulator transitions at much smaller values of U (~ Q /^jjsmgie-stte ^j^g^ 
density A^ = 1). However, the physics of the metal-insulator transition observed in the 
cluster-DMFT methods remains imperfectly understood, in part, because of the substantial 
computational expense of studying fully quantum cluster models. 

In this paper we present results of a systematic cluster DMFT study of a computationally 
simple polaron model, involving electrons coupled by an on-site density coupling to a disper- 
sionless classical oscillator. The model exhibits spectra and metal-insulator phase diagrams 
similar to those of the interacting electron models of primary physical interest, but allows 
an extensive analysis. We are able to determine the relation between short ranged correla- 
tions and the metal insulator phase boundary; elucidate the role of poles in the self energy, 
determine the minimum cluster size needed to see the effects of short ranged interaction and 
clarify the origin of the differences between dynamical cluster approximation (DCA) and 
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cellular-DMFT approaches. 

The rest of this paper is organized as follows. In section II we present the formalism 
we will use. We introduce the quantity "partial density of states" (PDOS) which is very 
important in determining the importance of short ranged effects. In section III wc review 
the polaron model, give some calculational details related to the impurity cluster problem, 
and discuss different theoretical ways to measure the strength of short-ranged correlation. 
In section IV we compare the results of one, two, and four-site cluster approaches applied 
to the half-filled model. In particular we show that the momentum dependent poles in self 
energies is the key to reduce the critical gap-opening U. Section V concerns how the poles in 
self energies depend on the temperature and interaction strength. In section VI we discuss 
doping dependence and give a criterion on the minimal cluster size required to capture the 
short-ranged charge ordering effects. In section VII we compared the results from DCA and 
cellular-DMFT, and from manganite and cuprate bands to illustrate the role of the partial 
density of states in cluster-DMFT calculations. Finally in a conclusion we summarize our 
results and discuss their implications. Two appendices give details of partial density of 
states and strong coupling calculations. 

II. GENERAL FORMALISM OF CLUSTER DMFT 

A. Overview 

In this section we review the general formalism for cluster dynamical mean field theory 
and give the equations for two frequently used approaches - cellular DMFT (CDMFT) and 
dynamical cluster approximation (DCA), and define a "partial density of states" (PDOS) 
used in subsequent analysis. Special attention is paid to the 2-site and 4-site DCA formalism 
on the square and cubic lattice. The main focus of this section is on the self-consistency 
condition. 

B. Cluster dynamical mean field formalism 

In an A^^-site cluster DMFT calculation, one solves a quantum impurity model which may 
be described as an A^^-site cluster coupled to a non-interacting electron bath. The solution 
is expressed in terms of the impurity Green's function and self energy - (5*"*^, E*"*p which 
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are Ns x A^^ matrices in the space of cluster labels. The dynamical mean field step requires 
relating components of G*"*^ and to some combinations of lattice Green's functions and 
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self energies. There are two frequently used algorithms: the cellular-DMFT approach [8|] 
which is defined in real space and the DCA approach [d^ defined in momentum space. 

In the cellular-DMFT to a given lattice, one chooses an A^^-site supercell which period- 
ically tiles the lattice and writes the Green's function as a matrix in this supercell, with 
entries depending on momentum k in the reduced Brillouin zone (RBZ) of the supercell lat- 
tice. One identifies sites of the quantum impurity model with sites of the supercell and take 
interactions of the quantum impurity model to be those of the original model, but confined 
to a supercell. The self consistency equation is 

G'^^p = [ {dk) L + /i - - t'M (1) 
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with [dk) = [d'^k) / {2txY ^ RBZ standing for the reduced Brillouin zone, and the A^^ x A^^ 
hopping matrix expressed in the basis of enlarged cell. In this approach the self-consistency 
equation is G-™'^ = G^'^*{Ri — Rj) with {R} basis vectors of the enlarged cell. The cellular 
DMFT method breaks the lattice translational symmetry because the intersite correlations 
(expressed in S*™'^) within a supercell are different from those between supercells. After the 
cellular-DMFT equations are solved one may restore the translation invariance by periodizing 
the self energy. Different methods have been proposed [l^ ; the issue is however not relevant 
to the considerations of this paper. 

The DCA approximation is constructed in momentum space. One tiles the Brillouin zone 
into Ns equal area non-overlapping tiles which we label by the average momentum Ki. One 
approximates the self energy S(A;, u;) by (cj) if k is in the tile corresponding to Ki. The 
impurity cluster model is defined by states labelled by the discrete set of labels Kf, the 
impurity model interaction is taken to be the momentum-space lattice interaction ^^^^^^^3^4 
but with wavevectors restricted to the set Ki. In the i^'j-basis the impurity model is diagonal 
and the self consistency equation is 

Gj'{u) = NsX [ (dp)— ^ ^ , , (2) 

Similar to single-site DMFT, it is useful to define the partial density of states (PDOS) for 
each Ki sector as 



A(e) = Ar, X / ^{dp) 5{e-e, 
JpeKi 



(3) 



(note that the prefactor A^^ leads to J deDi{e) = 1). Once the PDOS is given, the self 
consistency equation is 

G':^P{iu) = f de (4) 
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In all cases we have studied, the Ki basis diagonalizes the cellular DMFT impurity model 
(G and S) also. Therefore the cellular DMFT self consistency condition can be expressed in 
the form of Eq([2]) but the PDOS -Dj(e) will be different. An example is given in Appendix 
A. 

One notices that if all partial density of states are identical to the total density of states, 
the solution in single site DMFT is a solution in DCA, i.e. S"^'^ is identical for all Ki (S 
is purely local). A natural expectation is therefore that the separation between different 
sectors of PDOS indicates the degree of the short-ranged correlation. In the subsequent 
sections we show that this is indeed the case. 

While the DCA is naturally formulated in the cluster momentum basis Ki, it is often 
convenient to define a real-space basis by introducing a set of vectors {R}i=ir^Ns satisfying 

5,^ = ^f^e^i^^-^^y^' (5) 
■"^ 1=1 

(in this basis the on-site interaction term in the impurity model is identical to that in the 
lattice) Note that the {Ri} are only determined up to a uniform translation {Ri} — > {Ri+A} 
with A a lattice vector. Eq([5]) enables one to get real-space components from the impurity 
problem by the discrete Fourier transform 

Rm ]\[ — ^ 

i=l 
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GfP{uj) = -^\"e-'''^-'^-G'"'P{K,,u) (6) 



i=l 

where S*™P(i?m) = S*™;^^ and G^'^^^Rm) = G^^^^. For ith. momentum sector, the self 
energy is Tl'^'^{uj). 

In cellular DMFT, different choices of A^^-site supercell are possible. In DCA different 
tilings of the Brillouin zone are possible. A comparative survey of the different choices has 
not been undertaken and would be worthwhile. It is natural to expect that tilings should 
respect the point symmetry of the lattice as much as possible, i.e. if the lattice remains 



invariant under the transform r — > f = U-r, the ideal tihng of the Brillouin zone obeys 
E(fc) = T,{U^^k). In this paper we consider 2-site and 4-site tihng on square and cubic 
lattices plotted in Fig([T]). 

Notations: Since there will be many subscripts floating around in this paper, we explain 
the logic of our notation. For each quantity, like Green's function G, self energy S, and Weiss 
function a, there are momentum and real-space components labeled by some subscript. In 
this paper, the real-space component is labeled by a number (O:on-site, l:nearest neighbor..) 
while the momentum-space sectors labelled by capital letters (S, P, D..). 
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FIG. 1: (Color online) Partition of the Brillouin zone, (a) 2-site DCA on square lattice, (b) 4-site 
DCA on square lattice, (c) 2-site DCA on cubic lattice. 



C. Formalism for 2-site DCA 

Now we apply the machinery just described to specific cases, first to the 2-site DCA in 
the square lattice. After solving the 2-site impurity problem, one gets (in the real space 
basis) 
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(7) 



The partitioning of Brillouin zone in this case is given in Fig([T]), so two K points according 
to this division is Kj = 0, Kji = (7r,7r). We label region I and // or S and P sectors. 
Corresponding to Kj and Kjj, one gets Rq = and Ri = (±1,0) or (0, ±1) from the 
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requirement EqdS]). The lattice self energy is related to by 

T.DCA{k,Uj) = 

The partial density of states are 



= So + El for k G Region I{S) 
T!p^ = So - Si for A? G Region 11 {P) 



D^iP^{e) = 2 X 
and the self-consistency equation from Eq([2]) is 



s{p)\e) = zx I (dp) 5{e-eji) 
pei{ii) 



Go = (Gs + Gp)/2 
Gi = (Gs - Gp)/2 



with 



Gs{P) 



^5(P)(e) de 



uj + fx — ef. 



D. Formalism of the 4-site DCA Method 



(8) 



(9) 



(10) 



(11) 



In the 4-site DCA the Brillouin zone is divided into four sectors which are labelled as S, 
P, and D, as shown in Fig([T])(b). Four K points are (0, 0) (vr, 0) (0, vr) (vr, vr) leading to four 
R as (0, 0) (1, 0) (0, 1) (1, 1). The partial DOS is defined as 

ie)=Ax I (dp)6ie-e^) (12) 

Jp&S{P,D) 

where the superscript (4) is used to distinguish from the partial DOS in 2-site DCA. After 
solving a 4-site impurity cluster problem, in the disordered phase one gets 

^ Gq Gi G2 Gi \ ( So Si S2 Si ^ 
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G2 Gi Gq Gi 
\ Gi G2 Gi Go / 

and the momentum-dependent self energies are 
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(13) 



(14) 



and correspondingly the components of lattice Green's functions are 



Gs{P,D) 

The self-consistency equations are 



D 



S{P,D)\ 



Go = {Gs + 2Gp + Gd)/^ 

Gi = {Gs - Gb)/4 

G2 = {Gs-2Gp + Gn)/4: 



(15) 



(16) 
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FIG. 2: (Color online) The PDOS for 2-site and 4-site DCA partitioning on the square lattice with 
nearest neighbor hopping. The total bandwidth is 12. 



E. The band structure and partial density of states 

Throughout this paper we will consider the square and cubic lattices with only the nearest 
neighbor hopping and with total bandwidth chosen to be 12. For the square lattice with 
one S'-like orbital per site, the band dispersion is 

eg = -2t{cos{k^) + cos{ky)) (17) 

The corresponding PDOS for 2-site and 4-site DCA with bandwidth 12 (t = 1.5 in this 
case) are shown in Fig([2]). For the cubic lattice with one 5'-like orbital per site, the band 
dispersion is 

eg = —2t{cos{kx) + cos{ky) + cos{kz)) (18) 
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The manganite material llj are characterized by a pseudocubic lattice but because the im- 
portant orbitals are transition metal Cg orbitals for which the hopping is direction dependent. 
The band dispersion is described by a matrix as 



% = eoe + e^f^ + e^.f^ (19) 

where e, f^, and are unit and Pauli matrices; eo = —t{co?,{kx) + cos{ky) + cos(/c^)), 
= — t[cos(A;2) — {cos{kx) + cos{ky)) /2], and ex = —t^{cos{kx) — cos{ky)). The PDOS for 
2- site cubic lattice are plo' ' ' ' ^' Q 
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Frequency 

FIG. 3: (Color online) Partial DOS for 3D 25 and eg bands. 




III. THE POLARON PROBLEM AND ITS IMPURITY CLUSTER 



A. Overview 



In this section we briefly review the breathing mode polaron model and its solution 
by dynamical mean field theory. Based on the properties of this model, we provide two 
characterizations of the importance of short-ranged correlation. 

The breathing mode coupling: The Hamiltonian for the polaron model is 



Q2 



(20) 



where a, (5 can be spin or orbital indices, Qi is the local breathing mode distortion with U 
the restoring energy, and Ui = V ^, cj^Cia is the local charge density. This model has been 
solved using S-DMFT in Ref(l2|]. The main physical effect of the breathing mode Qi is 
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to trap electrons. In particular when the coupling U exceeds a density dependent critical 
Uc{n) the system becomes insulating at zero temperature. The physical interpretation of 
the insulating behavior is as follows: because a site with a higher density induces a lattice 
distortion which binds the electron and lowers the energy, the system eventually is phase 
separated into empty and maximally occupied sites where at empty sites there is no lattice 
distortion while at maximally occupied site the local distortion is at its strongest value. 
The insulating behavior is a consequence of the classical phonon approximation. Quantum 
fluctuations of phonons would lead to a small but non-zero conductivity at non-rational N . 



B. Solution by dynamical mean field theory 



Local partition function: In any cluster DMFT calculation, the impurity cluster problem 
is defined by a partition function. For the polaron model on an A^^ site cluster, the local 
distortion fields are described by Ng classical fields Qi which couple to a local density which 
is quadratic in fermionic operators. In this case, the fermionic degrees of freedom can be 
integrated out and the resulting partition function is an integral over Ng classical fields 



-V{{Q^})/T 



with 



V{{Qi}) = -2T X IVlog[a - diag(Qi, Q2, QatJ] + 3^7 E ^ 

i=l 



(21) 



(22) 



Here, a is a. Ns x Ns matrix composed of a set of Weiss functions {aj}. The 2 in front of the 
trace accounts for two degenerate local orbitals. 

For a 2-site cluster in the disordered phase, two sets of Weiss functions Oq, Oi are needed 
and 



i 2— site 



ai ao 

For a 4-site cluster, three sets of Weiss functions oq, oi, 02 are required and 

/ ao Oi a2 Qi ^ 
oi ao oi 02 



(23) 



^4:~site 
(Jj 



02 Oi Oq Oi 

y Qi 02 tti Qq ) 



(24) 
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The Weiss function, Green's function and self energies in real and imaginary axes are cal 



culated using the procedure described in Ref 
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131]. The main quantity we will show is 



the spectral function A{uj). The zero in frequency u is the Fermi energy. We define metal 
and insulator by the absence or presence of a gap at the Fermi level. To facilitate later 
comparisons, the total bandwidth is fixed at 12 for all band structure considered in this 
paper. 

Defining the importance of short-ranged effects: How to determine the degree of short 
ranged correlation is somehow arbitrary. Here we provide a definition in the context of 
cluster dynamical mean field theory: because the single-site DMFT completely neglects the 
spatial correlation, we define the strength of the short ranged effect from the deviation 
of cluster results from single-site DMFT results. Two quantitative measurements are now 
stated. First, both single site and cluster C-DMFT require a critical (minimum) interac- 
tion strength, f/i-''**^ and f/^'"«*e»~^ to get the insulating phase at zero temperature, and 
the ratio f/c'"«*er^f|i-s«te generally) is a measure of the importance of short ranged 

correlations. Second one can calculate the intersite phonon correlation function (QiQj) 
{= JZo'^dQk [QiQj] e-^(W^»/^/Z^™P). In single-site DMFT, {QiQ,) = {Qi){Qj) = while 
in cluster DMFT (QiQj) 7^ 0). Since that the number of intersite correlation functions 
increases exponentially when going to larger cluster size (including multi-site correlations), 
we shall primarily use the reduction of Uc in this paper. 



C. Pole in self energy 

We shall see that in the models we solve the insulating phases are associated with pole-like 
structures in the self energy 

S (25) 

uj- A-i'j ^ ^ 

characterized by amplitude V'^, position A, and damping 7 which is small or zero. The 
presence of a pole is related to insulating behavior because insulating behavior is character- 
ized by a gap, i.e. a region around the Fermi level (which we take to be 0) around which 
Im[G{k,uj)] = 0. If Im[G] = then /m[S] = except at poles. 

A gap in the spectral function at uj implies both that Im[T,{uj)] and that there is no 



11 



solution to Re[G = 0. In the DCA approximation this equation may be written as 

^-e^-^e[S^^M] = (26) 

for each momentum sector Ki and for in the allowed range corresponding to sector Ki. If 
the main contribution of the self energy is the pole, then the boundaries of the gap region 
can be obtained from the roots oi uo — t — V'^ / (cj — A) = which are 

^+(-) = i[(e + A) + HV(e-A)2 + 4F2] (27) 

Clearly a near Fermi level pole in E implies a large magnitude of -Re[E] at low uo, meaning 
that there is no solution to Eq (!26l) for the allowed range of e^. There are two interesting 
cases. If we consider density wave order at wavevector Q and neglect the long-ranged order 
by focusing on the diagonal elements of G, we find 

Sw(p,cu) (28) 

^ - V+Q 

(A broadened version of S ow is the basis of the Lee- Rice- Anderson theory of the pseudogap 
in thermally disordered charge density waves). 

S^jVK expresses the physics of level repulsion: if ^^j^q > then the presence of S_dvk 
pushes the solution of Eq (l2^ (which would he u = e^) to an energy lower than e^, while 
also creates a second solution at > ^^_^_q- We shall see that in many cases the cluster DMFT 
theories, which replace self energies by some sort of momentum averaged self energies, lead 
to just this physics, with the self energy corresponding to a i^" sector with dominantly 
less than having a pole at an energy A > 0, so that the effect of S is to push the states 
in this K sector away from ~ 0. This self energy arises from near resonant coupling to 
states in the complementary K sector. We interpret this as the density wave (Slater-like) 
mechanism. 

On the other hand, one sometimes finds cases where the self energy corresponding to a 
given K sector has a pole at an energy A in the middle of the allowed energy range arising 
from coupling to all K sectors. In that sector, the pole in self energy pushes states away in 
both directions. We identify this as the Mott mechanism. 

Extracting pole parameters: The pole parameters - the position, the amplitude, and the 
damping - from S ~ V^^/(ci; — A + 27) are obtained as follows. The position A is estimated 
by i?e[S] = (any Hartree contribution is absorbed into the chemical potential, which we 
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define to be 0); the amplitude V'^ and damping S are determined by fitting the i?e[E] to 

V^{iu-A)/[iu-AY + ^% 



IV. THE HALF-FILLED CASE 



Overview 



In this section the DCA method with cluster size 1, 2, 4 is used to analyze the two 
dimensional square lattice with nearest neighbor hopping t = 1.5. The model displays a 
metal insulator transition as the interaction U is increased beyond a critical value Uc, which 
as shown in Table I is in the range of 3 ~ 5 depending on the cluster size. 



cluster size 



1-site 



2-site 



4-site 



3.5 



Table I: Gap opening Uc for the polaron model on a square lattice with bandwidth 12 
obtained from DCA with different cluster size 

The excitation spectra are obtained for different cluster sizes. The physics of the metal- 
insulator transition and the variations in results with cluster size are related to the behavior 
of poles in the self energy. In the rest of this section we present and analyze the results for 
1-site, 2-site and 4-site DCA calculation for the local interaction range U < u'^-site ^jji-site 
is the gap opening U for single-site DMFT calculation). 

B. Single-site results 

We first present the single-site DMFT calculation. Fig(jll)(a) shows the spectral functions 
for U = 3,4,5 at T = 0.1. The gap-opening interaction for the single-site calculation, 
jji-site ^ is slightly less than 5. As discussed in the previous section, the insulating behavior 
is related to a pole in the self energy. This is shown in Fig(jl])(b) which reveals at t/ = 3 
the real part of the self energy is very weak; at t/ = 4 a strongly broadened pole begins to 
appear, and at [/ = 5 an undamped pole occurs. 

The self energy sum rule for this model is J ^/m[S(co')] = U"^. By fitting S (not shown 
here) at f/ = 5 we find the pole strength V"^ = 7.5 so that the pole contains roughly 30% 
of the total self energy sum rule. As U is increased the fraction of the self energy sum rule 
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Frequency Frequency 

(a) (b) 

FIG. 4: (Color online) Results from single-site DMFT calculation for two-dimensional square lattice 
with nearest neighbor hopping t = 1.5 (bandwidth 12), polaronic coupling (Eq (|20p ) f7 = 3, 4, 5 and 
T = 0.1: (a) spectral functions; (b) real part of self energies. 



contained in the pole increases to 90% of the total at f/ = 10 and 99% at f/ = 20. However 
even for U = 5 where there is a substantial non-pole contribution to Jm[E], the size of the 
gap is to a good approximation given by inserting E ~ 7.5 /u into the quasi-particle equation 
Eq([26]). 



2-site results 




Frequency 



(a) (b) 

FIG. 5: (Color online) Results from 2-site DCA calculation for two-dimensional square lattice with 
nearest neighbor hopping t = 1.5 (bandwidth 12), polaronic coupling (Eq (|20p ) {7 = 4 and T = 0.1. 
(a) Real part of self energies for S and P sectors, (b) Partial spectral functions for S and P sectors. 



In this subsection we show results from 2-site DCA calculation. We focus on U = 4 
where the solution is insulating but in the 1-site approximation would be metallic. Fig([5]) (a) 
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shows the real part of the self energy for 5* and P sectors at f/ = 4 T = 0.1 in for 2-site 
DCA method. The most prominent feature is poles. Comparison to Fig([2]) shows that the 
pole in the S sector (here A ~ 1.5) lies above the energy range where the S sector PDOS 
is nonvanishing (Fig([2])); while pole in the P sector lies below the support of the P sector 
PDOS. As discussed previously, the pole pushes the main part of the S sector PDOS down 
in energy while generating some excitation at energies above the pole energy. The opposite 
effect occurs in the P sector. Therefore in the 2-site DCA the origin of the gap opening 
transition is a coarse-grained (in momentum space) version of the Slater physics. 

D. 4-site results 




Frequency Frequency 

(a) (b) 

FIG. 6: (Color online) Results from 4-site DCA calculation for two-dimensional square lattice with 
nearest neighbor hopping t = 1.5 (bandwidth 12), polaronic coupling (Eq (j2U|) ) {7 = 4 and T = 0.1. 
(a) Real part of self energies for S P and D sectors, (b) Partial spectral functions for S P and D 
sectors. 

As in the previous subsection we focus on the insulating solution. In the 4-site DCA 
calculation, there are three different sectors S P and D in momentum space. Their self 
energies and partial spectral functions are shown in Fig([6]). As seen from Fig([2]), the PDOS 
of 5* and D sectors barely touch each other and the poles obtained from 4-site DCA have 
exactly the same effect as the S and P sectors in the 2-site calculation - S— and D— sectors 
repel each other. However in the P sector the pole energy is zero and the physics is "Mott- 
like" . We emphasize that in DCA equation all three sectors are coupled so states in the P 
sector are not isolated from states of other two sectors. Thus while the P sector bandwidth 
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is a factor of two smaller than the total bandwidth, the 4-site f/f--cc^ ^ 3.5 jg about 70% 
of the single site value. 

We observe that the polaron model (local coupling Qh+^Q^U, Eq ipl]) ) is in this respect 
different form the Hubbard model, where f/4-««te ^ OAU}''^^^'^ P, H]- We further observe that 
the [/-driven transition in the 4-site DCA approximation to the polaron model is continuous 
while in the Hubbard model the transition is strongly first order. In the Hubbard model 
the transition was associated with the dominance of a particular plaquette singlet state ^ 
not important here, suggesting that the smaller [/^-■^^te the Hubbard model case is a 
manifestation of the strong short-ranged ordering effect. 



E. Comparison between 2-and 4-site results 



T 




Frequency Frequency 

(a) (b) 

FIG. 7: (Color online) Spectral functions calculated by 2-site and 4-site DCA approximations 
for the two-dimensional square lattice with nearest neighbor hopping t = 1.5, polaronic coupling 
(Eq (pO]) ) [/ = 3, T = 0.2 (a) and C/ = 4, T = 0.1(b). 



Now we compare results for 2-site and 4-site DCA calculations in Fig(I7]). At U = 2 (not 
shown) the spectral functions are similar. At f/ = 3 (~ f/^"****^ but < f/^"****^) the spectra 
are different; but by f/ = 4 the spectra have again become very similar. As U continues 
to increase the difference between 2 and 4-site calculations remain small, but both remain 
different from the results of the single site calculation: in the multi-site case the gaps are 
larger and the gap edge more sharply defined. This is a consequence of short-ranged order 
(see also Table II, third row) and is demonstrated in Fig(l8]) in which the results obtained 
from single-site DMFT and 2-site DCA are shown for f/ = 5 and U = 10. One sees that the 
sharp band edge remains at large U. 
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FIG. 8: (Color online) Spectral functions from S-DMFT and 2-site DC A calculation for the two- 
dimensional square lattice with nearest neighbor hopping t = 1.5, polaronic coupling U = 5 (a) 
and U = 10 (b) , T = 0.1 and = 1. For U = 10 only the upper Hubbard band is shown. 



V. T AND U DEPENDENCE 



In this section we discuss the temperature and interaction dependence of the poles in the 
self energies. As previously discussed, in both single-site and cluster DMFT calculation, the 
insulating behavior can be traced back to a pole-like structure in self energy. The pole is 
characterized by its location A, amplitude V"^, and damping 6 which we now analyze for 
different interaction strengths U and temperature T. 

Fig® shows the temperature dependence of the density of states and the real part of 
the 5* sector self energies calculated in the 2-site DCA approximation for f/ = 4 (bandwidth 
12). As the temperature is raised, the gap begins to fill in, and this change is accompanied 
by an obvious decrease in the size and sharpness of the pole structure. While the pole is 
not well separated from the other contributions to S, a simple fit yields a gradual decrease 
in the amplitude (Eq fl2Sl) ) from 13 at T = 0.1 to 9 at T = 0.3, and a rapid increase in 
the damping 7 from negligible at T = 0.1 to 7 = 0.3 at T = 0.2 to 7 = 0.6 at T = 0.3. 
The shift in pole position is very small. Thus the main contribution to the destruction of 
the insulating phase is a broadening of the pole. A similar conclusion is found in the 4-site 
calculation (not shown here). The T-independence of the pole position is in contrast to 
the results of Refjl^ who found the pole position was temperature dependent, tracking the 
cluster spin correlation function. We believe this is a difference between the strong coupling 
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(a) (b) 

FIG. 9: (Color online) Results from 2-site DCA calculation for the two-dimensional square lattice 
with nearest neighbor hopping t = 1.5 (bandwidth 12), polaronic coupling U = 4 and T = 0.1-0.3. 
(a) Spectral functions, (b) Real part of self energies for S sectors. 



We next discuss the U dependence in DCA calculation. We begin with the strong coupling 
limit in which the excitations are expected to be roughly +{—)U. From Eq fl27l) . when 
amplitude of the pole is very large compared to both the pole location A and the half- 
bandwidth of non-interacting spectrum A, the excitations are roughly ~ +{—)V. 
Therefore in this limit one expects the pole amplitude is U"^ and the location of the pole 
is well within the gap, i.e. |A| << U (see the first row of Table II). As U is decreased, 
the bandgap and the pole energies both decrease, but the bandgap decreases much faster 
so that for the intermediate U the pole in the self energy lie actually within the continuum 
of excited states (see the first and second row of Table II), so is technically not a pole but 
is a resonance. Thus a pole in the self energy is not necessary for insulating behavior. We 
also find that the "pole" energy does not have a simple relation to the intersite phonon 
correlation {Q1Q2) (third row of Table II). This is again in contrast to results of Ref jl^. 
which found via an approximate analytical calculation that in the Hubbard model the pole 
energy was proportional to the spin correlation function {Si ■ 6*2). We believe the difference 
arises from the coupling strength. The U dependence is summarized in Table II whose first 
row gives the S sector pole position obtained from 2-site DCA calculation for different U at 
T = 0.1, second row the gap edge, third row the intersite correlation function. 
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u 


3 


4 


5 


10 


15 


S— pole location 


1.21 


1.31 


1.51 


2 


2 


gap edge 


~ 


1 


1.74 


6.75 


11 


-{QiQ2)/U^ 


0.40 


0.61 


0.73 


0.92 


0.924 



Table II: the pole position, gap edge, and intersite phonon correlation for different 
polaronic coupling U obtained from 2-site DCA calculation on the two dimensional square 

lattice with bandwidth 12. 



VI. FILLING AND CLUSTER SIZE 

A. Overview 

In this section we discuss band filling effects in cluster DMFT theory. The basic finding 
is that the filling determines the minimum cluster size required to provide a significant 
difference between single-site and cluster DMFT results - the minimum size increases as 
the system is doped from the half filling. In this section we restrict ourselves to the square 
lattice only. 

B. General filling 2-site 




Frequency 



(a) (b) 

FIG. 10: (Color online) Results from single-site and 2-site DCA calculations for the square lattice 
with nearest neighbor hopping t = 1.5, polaronic coupling U = 5, T = 0.1 at = 0.2. (a) Real 
part of self energies for S and P sectors, (b) Spectral functions. 
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In this subsection we compare the doping dependence found in single-site and 2-site 
DCA for U = 5 above the critical value for insulating behavior in both approximations. 
The discussion of Fig([S]) in Sec IV showed the pronounced differences between single-site 
and cluster DMFT calculations at density = 1. For the very low doping = 0.2 on the 
other hand, the single-site DMFT and 2-site DCA behave very similarly. Fig fllOp (a) shows 
that within DCA poles in self energies in the two sectors are almost identical, and are quite 
close to that in single-site DMFT. Consequently the spectral functions obtained from both 
methods are very similar ( Fig (fTOl) (b)). This can be understood as follows. At A^ = 0.2, the 
Fermi energy only lies deep inside the S'-sector PDOS, the P sector is irrelevant, and the 
DCA calculation within the S'-sector is just like a single-site DMFT calculation. 

C. Near quarter filling, 2-site and 4-site 



0.14 




Frequency Frequency 

(a) (b) 

FIG. 11: (Color online) Spectral functions calculated by single-site DMFT, 2-site and 4-site DCA 
approximations at [/ = 4, T = 0.1, N = 0.5(a) and N = 0.4(b) . 

In the classical polaronic model, insulating behavior can occur at any density if the 
interaction is strong enough. In this subsection we show that the Slater mechanism of shifting 
of states in k space can lead to insulating behavior at smaller U, if the band filling is such 
that one of the sectors is mostly filled, while the adjacent one with higher energies is nearly 
empty. However these cluster size effects diminish rapidly in importance as the density is 
moved away from commensurate values. The left panel of Fig ffTTl) compares the results from 
1-site, 2-site, 4-site DCA approximations to the two dimensional square lattice with nearest 
neighbor hopping at density A^ = 0.5 and intermediate U. The four site approximation 
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(a) (b) 

FIG. 12: (Color online) Self energies calculated by 2-site (a) and 4-site (b) DCA approximations 
for the two dimensional square lattice with bandwidth 12 at C/ = 4, T = 0.1, = 0.5. 



produces a clearly insulating behavior, the two site gives a moderate "pseudogap" and the 
single site approximation yields results very similar to the non-interacting model. The right 
panel of Fig fllll) shows that the differences decreases rapidly as the doping is moved from 
the commensurate value. 

To unravel the origin of the insulating behavior we consider the self energies for the 4-site 
DCA shown in Fig (fT2|) . The right panel shows a clear pole in the S sector, lying at an 
energy slightly above the chemical potential fi and substantially above the average energy 
of 5* sector fi — 0.93. This pole gives rise to the Slater effect of pushing most of 5* sector 
states down in energy. Consequently the P and D sector self energies lie below the P and 
D PDOS and push these states up, leaving the gap observed in Fig (fTTi) . 

By contrast the left panel of Fig (fT2|) shows for the 2-site DCA an interesting hybrid 
behavior. The P sector pole lies below the P-sector PDOS and pushes these states up. The 
S sector pole lies roughly at /i — 1, below the average energy of S PDOS (e)!"''**'^ = /i + 0.95 
and acts in a Mott-like fashion, splitting the band into two. 

VII. THE ROLE OF PARTIAL DENSITY OF STATE 

In this section we establish the relation between the PDOS defined in Eq(l3]) and the 
differences between cluster and single-site DMFT calculations. In essence, the more overlap 
between PDOS, the more similar are the single-site and cluster calculations. To measure the 
overlap between PDOS of two sectors, we first compute the overlap Oij = J deDi{e)Dj{e). 
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One can also use the first two moments of tlie partial DOS. For PDOS i, (e)j = J deeDi{e) 
gives the average position of the partial DOS while ai = a/ (e^)i — (e)f tells how spread the 
PDOS is. The overlap is related to {e)i/ai. We consider two examples. First we compare 
2-site DCA and 2-site cellular DMFT solutions to the polaron model on the square lattice. 
Second, we compare the approximations to two models involving the polaron model defined 
on a cubic lattice: one has two degenerate S'-like orbitals per site the other two Cg orbitals, 
with dispersions given in Eq ffTSl) and Eq f|T9|) respectively. 



A. Difference between 2-site DCA and 2-site CDMFT for 2D square lattice 
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FIG. 13: (Color online) Partial DOS for 2-site DCA and Cellular-DMFT on the square lattice with 
nearest neighbor hopping (t = 1.5 (bandwidth 12). Uc for 1-site DMFT, 2-site cellular DMFT, 
2-site DCA approximations are roughly 5, 4, 3 respectively. 
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In this subsection we compare the two site DCA and cellular DMFT solutions to the 
polaron model defined on the square lattice. As discussed in Sec II and Appendix A the 
mean field equations may be expressed in the same form as 

G™^ = [ dt ^^'^j'^ , , (29) 

The two approximations are distinguished by the differences in PDOS, shown in Fig( fT3ll . 
In the DCA case two PDOS do not overlap (the small apparent overlap is the result of 
numerical broadening); in the cellular DMFT the overlap is seen to be more substantial. 
Quantitative information on the overlaps is given in Table III 
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{e)s = - 




(Ts=ap 


DCA 





-2.41 


1.82 


Cellular-DMFT 


0.15 


-1.48 


1.89 
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Table III: The overlap information of PDOS of 5* and P sectors for DCA and cellular 
DMFT on the square lattice with bandwidth 12. {e)s{p) and as=o'p are the average energy 

and the standard deviation of S{P) PDOS. 

To further explore the difference we show in Fig ffT^ the spectral functions calculated at 
U = 4 and at the relatively high T = 0.2. The DCA calculation reveals a well defined gap, 
whereas the cellular DMFT reveals a weak pseudogap which will evolve to a small true gap 
as T 0. As seen from Eg (1291) these differences can arise mathematically only from the 
difference in PDOS. 
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FIG. 14: (Color online) 2-site cellular-DMFT and DCA spectral functions calculated for the square 
lattice with nearest neighbor hopping t = 1.5 at half-filling, polaronic coupling U = 4, T = 0.2. 
The DCA result is more insulating (having smaller density of state around the Fermi energy) than 
cellular-DMFT. 

Solving the equations reveals that the gap-opening U of single-site DMFT at half filling 
is U^-^^te _ 5 (^^Q^g^i bandwidth is 12), for 2-site cellular DMFT ^/s-^^e^nZar _ 4 ^^^^ 2-site 
DCA, jj^-D'^^ 3. The difference shows that the DCA expresses short ranged correlations 
more strongly than the cellular DMFT. 

B. 3D eg and 25" band 

In this subsection we compare the results of applying the DCA approximation to two 
orbitally degenerate models defined on a cubic lattice: the S-orbital model, with two de- 
generate bands of dispersion given in Eq (fT8|) and the eg model, appropriate to the colossal 
magnetoresistance manganites, with dispersion given by the eigenvalues of Eq ([T^ . We 
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T 



SDMFT 

-— 2-site CDMFT 
I - - 2-site DCA 




choose the bandwidth to be 12 for both bands. The corresponding PDOS are shown in 
Fig ([3]) and the overlaps are given in Table IV. 
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{e)s (=-( 




(^s (= o-p) 


25*— like band 


0.044 


-1.84 


1.8 


Cg band 


0.052 


-1.82 


2.98 



Table IV: The overlap information of PDOS of 5* and P sectors for S (Eq ffTSl) ) and Cg 
(Eq (fT9|) ) bands on the cubic lattice with bandwidth 12. {e)s{p) and as=crp are the average 
energy and the standard deviation of S{P) PDOS. 

We see that the eg band has larger PDOS overlap than the S band and thus one expect 
a smaller short ranged correlation for the Cg case. Indeed, for 2S band (f/^"''**^ ~ 4) 
jj2-DCA^jji-site _ _ 0.75 while for Cg (t/]"^^*^ < 6) f/^-^^^/f/J - site ~ 5/6 = 0.83. 
The corresponding spectral functions are provided in Fig( fT5l) . We observed that in the 
insulating region the DCA leads to a larger gap and sharper band edge than the single-site 

DIVIFT but ^"^'^ rl iflFoTcin r^oc qto locc in "i^lnci ^ naca 
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(a) (b) 

FIG. 15: (Color online) Spectral functions for the coupling close to the gap opening Uc- (a) For 
two degenerate 5-band with bandwidth 12, the ~ 4 (solid) and U^~^'^'^ ~ 3 (dashed), (b) 

For Cg band with bandwidth 12, the Ul~^^^^ ~ 6 (solid) and U^~^'~^^ 5 (dashed). Those curves 
are calculated at T = 0.1. 

VIII. SUMMARY AND DISCUSSION 

We first summarize our findings. The cluster-DMFT calculation reduces the critical inter- 
action strength to open a gap. The magnitude of this is introduced as a quantitative measure 
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of the importance of short-ranged correlation. The partial density of states (PDOS) and the 
filling of the system are two crucial quantities in the cluster DMFT calculations. The reduc- 
tion of critical U in cluster calculations is mainly caused by the momentum dependent pole 
structure in self energies which pushes PDOS around Fermi energy away from each other. 
Similarly we find that significant differences between single-site and cluster calculations only 
occur when the Fermi energy involves several PDOS (2 at least). For example in the square 
lattice, the 2-site produces small U insulating behavior for half-fiUing problem but differs 
little from single-site DMFT for quarter-fiUing. Obtaining a small U insulating state in the 
quarter filling requires a 4-site cluster. We show that when the cluster size is large enough, 
a large separation between PDOS around Fermi energy implies a strong short-ranged effect. 
One notes that the difference between single-site and cluster DMFT results becomes pro- 
nounced for the range of U where < U^~^^^^, above which the cluster calculations lead to 
a larger gap and sharper bandedge. These features persist in the strong U limit and are 
consequences of the short-ranged correlation included in the cluster approximation. 

The use of cluster DMFT methods to study metal insulator transitions at carrier con- 
centrations different from half filling is an intriguing issue. Our results suggest that the 
cluster size must be tuned according to the filling to be studied. Roughly filling l/Ng re- 
quires an A^5-site cluster. Thus for example study of the interesting phenomena associated 
with "half-doped" manganites (one carrier for every two two-fold degenerate Mn orbitals, 1 
electron per 4 orbitals, quarter filling), at least a 4-site cluster is required. The physics of 
the metal insulator transitions at 7^ 1 is found generally to be Slater-hke, arising from 
the separation in energy of k space sectors. 

A generic results is that A^^-site [Ng > 1) cluster DMFT calculations predict smaller 
critical U for the gap opening (metal-insulator) transition than does the single site approx- 
imation. Although various effects occur, we find that generically the most important role is 
played by a Slater-like effect. The calculation is most naturally viewed in terms of momen- 
tum space sectors, and the dominant contribution to the insulator behavior is given by a 
pole structures in the self energy which act to open gaps between the different momentum 
sectors. 
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APPENDIX A: PDOS OF CELLULAR DMFT APPROXIMATION 



In this appendix we derive the partial density of states for 2-site cellular DMFT approx- 
imation for square lattice with nearest neighbor hopping t. We begin by casting these two 
methods in the same formalism. Once the 2-site impurity cluster problem is solved, one 
gets EqQ. The self consistency equation for both 2-site DCA and cellular-DMFT can be 
expressed in the same form similar to Eq([T]) as 




2 X / (dk) 
'kei 



id + n- 




(Al) 



Note that Region / is also the reduced Brillouin zone as depicted in Fig([T]) (b). Choice of 
ipj; distinguishes DCA and CDMFT. For DCA, ip^: = er 



= -2t(cos k^ + cos K); for CDMFT, 
tfj,: = (f),: = -t[l + e-"^'^- + e-'^'^^^-^y^ + e-^'^^^+^y^] (see Eq(lO) in RefpJ]). The partial 
density of states in this formalism is defined as 
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S{P) 



,e = - 



X / {dk) 





n 



-1 -1 



- 11(22) 



(A2) 



Note that the off-diagonal terms are zero for both cases. For ip^ = eg one retains the 



definition in DCA; for ipv 



one obtains the PDOS in cellular DMFT as 



Die)s(p) 



2 X / (dk) 6{e 

'k&I 



(A3) 



APPENDIX B: STRONG COUPLING, ZERO TEMPERATURE LIMIT IN 2- 
SITE DCA 



In this appendix we analyze the polaron problem in the large U , zero temperature limit. 
To facilitate the calculation, we assume that the partial density of states for S and P sectors 
to be semicircular with the average position at negative and positive a, i.e. 
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S{P){ 



a 



fBll 
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because of the analytical expressions. We will show that in this hmit, the self energies in S 
and P sectors indeed have pole as 

t/2 



(B2) 



u;-(+)A 
with the pole A — a. 

We first solve the impurity cluster problem. Following the notation in Sec.III, in the 
T = 0, large U limit, {Qi,Q2) are either (+[/,—[/) or (—[/,+[/). Averaging these two 
configurations, one obtains 



(B3) 



and interpreting 



1 



SiP) a^^{-)a,-U^/[a^-{+)a,\ 
The self energy is obtained by E = a — ((J*™?)-!. After a httle algebra, one gets 

E - 

^S{P) - 7-T 

ao - (+)ai 

Using this self energies, the lattice Green's function for S sector is 
Gfiuj) -- 



(B4) 



(B5) 
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Dsje) de 
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l/{u; + a - E5) 



(B6) 



Similarly Gp^{u) ~ l/{uj — a — Sp). Comparing the lattice and local Green's function, one 
gets ao + (— )ai — co + (— )q; thus 

C/2 C/2 



Oo - (+)ai <^ - (+)« 



(B7) 
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